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Abstract 

We calculate the antiferromagnetic spin wave dispersion in the half-filled (electronic 
density n = 1) Hubbard model for a two-dimensional square lattice, using the random 
phase approximation (RPA) in a broken symmetry (spin density wave) ground state. Our 
results for the spin wave dispersion, <jj(q), are compared with high-resolution inelastic 
neutron scattering performed on La2Cu04. The effects of different band structures and 
different values of the on-site Coulomb interaction on the spin wave spectrum is studied. 
Particular attention is put on the high energy dispersion values <j(7t/2, 7t/2) and uj(0, it). 

Introduction: In two recent papers, [1, 2] high-resolution inelastic neutron scattering mea- 
surements have been performed on two different two-dimensional spin 1/2 quantum antifer- 
romagnets. These are copper deuteroformate tetradeuterate (CFTD) and La2Cu04. Sur- 
prisingly, the dispersion at the zone boundary that has been observed in the two materials, 
does not agree with spin- wave theory predictions [3]. Moreover the amount of dispersion is 
not the same for both materials. In CFTD the dispersion is about 6% from lo(ti/2, n/2) to 
ui(n, 0), whereas in La2Cu04 it is about -13% along the same direction. In the case of CFTD 
the dispersion at the zone boundary can be explained using the nearest-neighbor Heisenberg 
model alone, [2] and high precision quantum Monte Carlo simulations have confirmed that 
it is so. [4] On the other hand, an explanation for the observed dispersion in La2CuC>4 has 
been proposed [1] using an extended Heisenberg model [5, 6] involving first-, second-, and 
third-nearest-neighbor interactions as well as interactions among four spins. 

In a previous paper [7], we have shown that it possible to obtain the observed dispersion 
difference of -13% for La2Cu04 using the single band Hubbard model at half filling, with 
nearest neighbor hopping. In our formulation the extended Heisenberg model used in ref. [1] 
is incorporated by means of the virtual excursions of the electrons on the lattice. Fitting our 
results to the experimental data the obtained values of U and t agree well with those of ref. 
[1] and where confirmed by Quantum Monte Carlo calculations in the Hubbard model [8]. 

In this paper we generalize our previous study incorporating in the calculations the effect 
of a second nearest neighbor hopping t' in the electronic spectrum (in high-T c materials the 
ratio \t' /t\ ranges roughly from 0.1 to 0.5). The effect of U on the spin wave dispersion at 
the special points uj(tt/2, tt/2) and Lj(0,n) is also studied. 
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Model Hamiltonian: The Hubbard model for a square lattice of N sites is defined as 

h = Yl^~^H,a c ^ +u J2 c t^ c ti Ci - 1 ' (1) 

h,a i 

where e(k) defines the energy dispersion for independent electrons. In this work we consider 
two different electronic energy dispersions given by 

e(k) = —2t cos k x — 2t cos k y (2) 

and by 

e(fc) = — 2tcosk x — 2t cos k y — At' cos(k x )cos(k y ) . (3) 

The first energy dispersion has the nesting vector Q — (tt, it); the second one is not nested. 

The broken symmetry state is introduced by considering the existence of an off-diagonal 
Green's function given by 

F a {p-r - t') = -< T T c ?±(3jCT (r)4 jCT (r') > . (4) 
in addition to the usual Green's function: 

G CT (p; T - t') = - < T T c ff Ar)cijr') > . (5) 
At the mean field level the Fourier transform of these two Green's functions are given by 



itu n - E + (p) itu n - (p) 
F„(p,iuj n ) = . ^ + . % R , (7) 

where the energies E± are given by 



E±(P) = m + ^ + Q) + V\ ± \sfm - a? + Q)? + IPm* , (8) 
= e(p) — fx, and the coherence factors read 

E + - gp + Q) - Un/2 E+-qp)~Un/2 

U P — 771 TTi > W P — Tp TP ' W 



and 



Uma/2 _ Uma/2 

V V,° = 1? F~ • ( 10 ) 



p ' a TP TP ' P ' CT _ E 1 

£/_(- — rj — £L/-^ 

Spin susceptibility and spin waves: In order to describe the spin dynamics of the system 
we consider the transverse spin susceptibility \ h(<7>* a 'n), which is defined as 



X -+{q,^n)=A / dre— T <T T S-(q,r)S + (q,0)>, 
Jo 



(11) 
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where (3 = 1/T is the inverse temperature, T T is the chronological order operator (in imaginary 
time), S~(q) = X)p c ]? | c p+9.T an ^ S + (q) = [S~ (q)]^ . The above expression can be written as 

X+-(q,iw n ) = H 2 b^2 dr^e lUnT 

r 

<T T [- / df J ff c/ (f)] n 4 i (T)c ?+9 - T (r)4, + . T (0)v, i (0)> d . c ., (12) 
Jo 

where d.c. stands for differently connected diagrams. The susceptibility is evaluated at the 
RPA level [7] and the spin wave spectrum Lo(q) is determined from the poles of x+-(<7, u>+iO + ). 
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Figure 1: Spin wave dispersion, in meV, along high symmetry directions in the Brillouin 
zone. The circles are the data reported in Ref. [1] at 10 K. The dashed line is the 
analytical result (at K) for U = 1.8 eV, t = 0.295 eV. The solid line is the analytical 
result (at K) U = 2.1 eV, t = 0.34 eV, and t'/t = -0.25. The momentum is in units of 
2tt and M = (1/2, 1/2), X = (1/2, 0), and T = (0, 0). 



In Figure 1 the spin wave spectrum is plotted along high symmetry directions of the 
Brillouin zone for both dispersions (2) and (3). For the case of the dispersion (2) the values 
of t and U giving the best fit to the experimental data are, for T = K, U = 1.8 eV, 
t — 0.295 eV, with U/t = 6.1 (we note that there is no measurable change between the 
calculation at T — K and T — 10 K). For the dispersion (3), which does not present 
nesting, the experimental data can be fitted using U = 2.1 eV, t = 0.34 eV, with U/t = 6.2, 
and t'/t = —0.25. These last set of values agree with those determined in ref. [9], from an 

perturbative calculation of the poles of Xh (<f> w + *0 + ). It is clear from these results that 

the introduction of a more realistic band structure, given by (3), leads to larger values of t 
and U. 
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Figure 2: Variation of ui(0, 7r) and ui(tt/2, it/2) as function of U/t. The lines with circles 
and squares refer to the case where t' /t = —0.25; the other lines refer to the case t' = 0. 

Let us now consider the effect of the Coulomb interaction U on the energy difference 
w(0, it) — u(-7r/2, 7r/2). Experimentally the dispersion, in CFTD, is about 6% from u)(tt /2,ir /2) 
to u)(tt,0), whereas in La2Cu04 it is about -13% along the same direction. The question is 
whether it is possible to obtain these very different behaviors from a single model Hamiltonian. 
In Figure 2 we study the evolution of uj(Q,-k) and uj(tt/2,tt/2) as function of U for the 
dispersions (2) and (3). For the dispersion (2) it is clear that lo(0, it) is always larger then 
lo(tt/2,tt/2), except in the limit U — > oo where they become equal. This behavior was 
confirmed by Quantum Monte Carlo [8]. On the other hand, for the dispersion (3), u>(0, tt) 
and w(7r/2, it/2) become equal for a finite value of U/t, which for the chosen parameters 
is U/t ~ 10.5. For larger values of U/t the dispersion at uj(tt/2,it/2) becomes larger than 

W(0,7T). 

In conclusion, our results shown that it is possible to fit the spin wave spectrum with the 
t — t' — U Hubbard model and that the two different behaviors for ui(0, 7r) and lo(tt/2,it/2) 
observed in La2Cu04 and in Cu(DCOO)2-4D20 follow from a single model Hamiltonian. 
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